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ABSTRACT 

Identification of high redshift clusters is important for studies of cosmology and cluster evolution. Using 
photometric redshifts of galaxies, we identify 63 1 clusters from the Canada-France-Hawaii Telescope (CFHT) 
Wide field, 202 clusters from the CHFT Deep field, 187 clusters from the Cosmic Evolution Sui-vey (COSMOS) 
and 737 clusters from the Spitzer Wide-area InfraRed Extragalactic survey (SWIRE) field. The redshifts of 
these clusters are in the range of 0.1 < z < 1.6. Merging these cluster samples gives 1644 clusters in the four 
survey fields, of which 1088 are newly identified and more than half are from the large SWIRE field. Among 
228 clusters of z > 1, 191 clusters are newly identified, and most of them from the SWIRE field. With this 
large sample of high redshift clusters, we study the color evolution of the brightest cluster galaxies (BCGs). 
The colors r' — z' and r+ — ma.g^m of the BCGs are consistent with a stellar population synthesis model in 
which the BCGs are formed at redshift Zf > 2 and evolved passively. The colors g' — z' and B — ms g^m of 
the BCGs at redshifts z > 0.8 are systematically bluer than the passive evolution model for galaxy formed at 
~ 2, indicating star formation in high redshift BCGs. 

Subject headings: galaxies: clusters: general — galaxies: elliptical and lenticular, cD — galaxies: evolution 



1. EVITRODUCTION 

Galaxy clusters are important objects to study the large 
scale structure (Bahcall 1988; Postman et al. 1992) and con- 
straint the cosmological parameters, e.g., fl„i, the mass den- 
sity parameter of the universe, and as, t he amphtude of mass 
fluctuations at a scale of 8 /t~^ Mpc dBahcall et alj 119971: 
iReiprich & Bohringer| |200l IWen et al.l 120101 ). Clusters are 
also important laboratories to investigate the properties of 
dark matter, hot gas, gal axies and active galactic nucleus in 
dense environment (e.g.,|Castillo-Morales & Schindlerll2003l : 
lDressledlT980t ICroft et al.ll2007h . However, most of these in- 
vestigations are based on clusters in the local universe. 

High redshift clusters can provide information on cosmo- 
logical structure in the vast universe and the evolution of 
cluster properties with cosmic time. It has been shown that 
cluster galaxies at high re dshift have more star fo r mation 
than those at low redshift ("B utcher & Oemledll978l |1984|) . 
The evolution of cluster mass functions gives constraints on 
not only the cosmological parameter s but also dark energy 
equation of state parameter loq (e.g., IVikhlinin et al.ll2009l) . 
The abundance of (even a few) high redshift massive clus- 
ters constrains the n on-Gaussianity of primordial perturba- 
tion dMatarrese et al . 2000; Hoyle et al. 2010). High red- 
shift clusters are needed to stu dy the form ation and ev olution 
of brightest cluster gal axies (iMiley et al.i,2006: .Whiley et alJ 
l2008tlStott etaPIIOOSh . 

Li. Previous detection of high redshift clusters 

Tens of thousands of galaxy clusters have bee n identified 
from various surveys i n the last decades (e.g., lAbell et aP 
[1981 IWen etaUlIOOl EaoetaDfeOlO). They mostly have 
redshifts z < 0.5. High redshift clusters were recently 
identified from multi-color optical or infrared deep surveys. 
Only a few hundreds of clusters have r edshifts z > 1 (e.g., 
iGladders & Yedl2005tlGotoetani2008|). 

Using cluster red sequence method, iGladders &Yeel(l2005l) 
detected 429 candidate clusters or groups with redshifts of 
0.2 < z < 1.4 from the Red-Sequence Cluster Survey (RCS), 
of which 67 have 0.9 < z < 1.4. Using the Cut and En- 



hance method, [Goto et all ( l2008l) identified 16 cluster candi- 
dates atredshifts_0^ < z < 1.7 from the AKARI deep sur- 
vey. lEisenhardt et al.l (|2008|) used a wavelet algorithm based 
on photometric redshift to the Spitzer Infrared Array Camera 
(IRAC) Shallow Survey and NOAO Deep Wide-Field Survey 
data and identified 335 cluster and group candidates, of which 
1 06 have reds h ifts z > 1. 

lOlsen et al] (|2007|) applied a matched-filter cluster de- 
tection algorithm to the Canada-France-Hawaii Telescope 
Legacy Survey (CFHTLS). They identified 162 cluster can- 
didates over an area of 3.1 12 deg^, of which 29 have redshifts 
z > 1 . iGrove et al.l (l2009l) found 114 cluster candidates by us- 
ing a matched filter detection method to the CFHTLS r' band 
data of the Deep fields, and 247 cluster candidates from the 
z' band data. Merging the samples in both bands gives 233 
cluster s, of which 93 have redshifts z > 1. iThanjavur et al.l 
(120091) found 5804 cluster candidates from the CFHT Wide 
field of 161 deg^ using the red sequence method, of which 
13 have redshifts z > 1. By mapping densit y of galaxy dis- 
tribution in photometric redshift space, Ada mi et aP (1201 Oh 
detected 1200 candidate clusters with masses greater than 
1.0 X 10" Mq from the CFHTLS Deep and Wide fields, 
of which 302 have redshifts z > 1. iMilkeraitis et al.1 (1201 Oh 
applied a 3 -dimensional matched-filter cluster detection algo- 
rithm to the CFHTLS data and identified 673 cluster candi- 
d ates in the red sh ift ra nge of 0.2 < z < 1.0. 

iKnobe l et al.l ((2009) used the friends-of-friends and 
Voronoi tessellation methods and identified 102 groups of 
z < 1 with more than five member galaxies from 10,000 
redshif ts from the Cosmic Evolu tion Survey (COSMOS) 
field. Ivan Breukelen et al] (|2006|) identified 13 clusters in 
the redshift range of 0.61 < z < 1.39 from the UKIRT 
Infrared Deep Sky S urvey (UKIDSS) Eariy Data Release. 
IZatloukal et alJ (l2007h identified 12 cluster candidates in the 
redshift range of 1.23 < z < 1.55 from the Heidelberg In- 
fraRed/Optical Cluster Survey (HIROCS) survey in the COS- 
MOS field. Chiaberge et al. (2010) found three candidate 
clusters of galaxies at red shifts most likely between 1 .7 and 
2.0 in the COSMOS field. IWang et all ([20 IQ) used the friend- 
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of-friend method to the Great Observatories Origins Deep 
Survey (GOODS) data and identified 206 group at redshift be- 
tween 0.4 and 1.0. Fro m the first 36 XMM-New ton pointings 
on the COSMOS field. iFinoguenov et alj ( l2007h identified 72 
clusters, of which 8 have redsh ifts z > 1. From the XMM- 
LSS survey. iPierre et all (l2007|) identified 73 clusters. From 
the Subaru-XMM Deep Field. Finogue nov et"an (|2010) iden- 
tified 57 cluster candidates, of which 13 have redshifts z > 1. 

Besides the sample finding, some authors have spectroscop- 
ically found or c onfirmed individual high redshift clusters. 
iKurk et al] (l2009l) confirmed a galaxy cluster of z = 1.6 
from the Galaxy Mass Assembly ultra-deep Spectroscopic 
Survey (GMASS). ^ilbank et al. (2008) confirmed a com- 
pact supercluster structure c omprising three clusters at 2 = 
0.9. iPapovich et alj (1201 Ol) discovered a galaxy cluster at 
z = 1.62 from the Spitzer Wide- Are a Infrared Extr a galac- 
tic survey (SWIRE) XMM-LSS field. lAndreon et all (120081) 
confirmed a cluster at z — 1.016 using a modified red se- 
quence method, follow up spectroscopy and X-ray imaging. 
From the Spit zer/IRAC Shallow Survey of the Bootes field, 
[Stanford et al.l (12005) confirmed a galaxy cluster z — 1.41. 
From the Spitzer Adaptation of the Red-sequence Cluster Sur- 
vey (SpARC S), Muzzin et al. (2009), Wilson et al. ( 2009) and 
lDemarco"et al.l (12010 ) confirmed two clusters at z ^ 1.2, 
one cluster at z = 1.34 and three cluster at z — 0.87, 
1.16, 1.21, resp e ctivel y. Using a 3-dimensional technique, 
ICastellano et al.l (1200 7*) detected a forming galaxy cluster at 
redshift 1 .6 in the GOODS field. A few high redshift mas- 
sive clusters were found by X-ray observations and Sunyaev- 
Zeldovich (SZ) effect. From the XMM-Newton observa- 
tions, a massive X- ray cluster at z — 1.39 was discovered by 
iMulHs et al.' ('2005') and a massive X-ray cluster at z = 1.45 
by Stanford et al. (2006), and a cluster at z = 1.22 was iden- 
tified by Bremer et al. (2006) and a cluster at z = 0.95 by 
iFassbender et a l. (2008). Using the SZ effect. iBrodwin et"al] 
(1201 Oh found a massive cluster at z = 1.07 from the South 
Pole Telescope (SPT) data. 

1.2. Evolution of brightest cluster galaxies 

The brightest cluster galaxy (BCG) is an elliptical galaxy 
located at the potential center of a galaxy cluster BCGs are 
the most luminous galaxies in the universe and in general 
have no prominent ongoing star formation. They are red in 
the rest-frame color Their surface brightness profiles are dif- 
ferent from those of ordinary elliptical galaxies (non-BCGs), 
and they do n ot follow the basic scaling relations of normal el- 
lipticals (e.g..' Lauer et al.ll2007l: lUiu et al.ll2008h . Because of 
the dominant role inside clusters and their unusual properties, 
the formation and evolution of BCGs are very intriguing. 

The properties and evolution of BCGs were studied using 
identified high redshift clusters. iBernardil (l2009l) noticed that 
BCGs at lower redshift have larger sizes and smaller velocity 
dispersions, indicat ing that early -type BCGs grow from many 
dry mino r mergers (pVIilev et al.ir200 6 ) rather than a few major 
mergers (lRinesetal.1120071) . By comparing the structural pa- 
rameters of BCGs of nearby clusters (0.04 < z < 0.07) and 
those at intermediate redshift (0.3 < z < 0.6). lAscaso et al.l 
(I2OII) confirmed the size decrease of BCGs with redshift. 
The properties of BCGs are related to the host clusters. More 
luminous, larger and more centrally located BCGs are located 
in more massive and rich galaxy clusters. 

On the -fC-band Hubble diagram, BCGs do not exhibit 
any luminosity evolution with r edshift (Collins & Mann 19981 
lAragon-Salamanca et al.|[T998t iWhilev et al.1 120081) and sig- 



nificant change in the stellar mass (IWhilev et al.ll2008h . The 
colors of BCGs are in good agreement with the evolved old 
stellar population formed at z > 2. This means that the stellar 

J iopulation in BC Gs has been in place since at lease z = 2 
Stott etal.l[2008h and the avera ge stellar rn ass of BCGs re- 
mains constant since z ~ 1.5. IStott et aH {2010) therefore 
concluded that dry merger seems to have little effect on evo- 
lutions of B CG stellar mass o ver the last 9-10 Gyr 

However, iLiu et al.l (l2009l) found that dry merger plays 
an important role in the stellar mass assembly of BCGs. 
iDe Lucia & Blaizo3 (|2007|) concluded from hierarchical sim- 
ulations that the half mass of a typical BCG is assembled by 
dry mergers after z ^ 0.5. The stellar populations of BCGs 
are formed very early, 50% at redshift z ~ 5 and 80% at 
z - 3. 

Observationally, a very small number of BCGs, e.g., those 
in Abell 1835, Zw3146 and MACS J0913.7H-4056, excep- 
tionally show the features of star formation (McNamara et al 
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starburst features (iLiu et al.lliol 1 ). The star formation is the 
dominant power for the infrared and Ha emission detected 
from BCGs (O'Dea et al. 2008). 

As shown above, many controversies exi st on the forma- 
tion of BCGs a nd the role of dry merger (Liuet al.ll2009[ 
IStott et al]l2010l) . It is not clear when the BCGs are formed, 
what process dominates the evolution of BCGs. Finding more 
high redshift clusters and their BCGs is crucial to study these 
questions. 

In this paper, we identify high redshift clusters from the 
CFHT Wide field, the CHFT Deep field, the COSMOS field 
and the Spitzer SWIRE field. In Section 2, we first apply our 
cluster detection algorithm to the photometric redshift data of 
the four fields, and identify high-redshift clusters. We then 
discuss the false detection rate, the accuracy of determined 
cluster redshift, cluster richness and X-ray luminosity of some 
clusters. The color evolution of BCGs is studied in Section 3. 
Conclusions of this paper are presented in Section 4. 

Throughout this paper, we assume a ACDM cosmology by 



taking Hq =100 km s 
and JIa — 0.7. 



Mpc~ , with h 0.72, r2„ 
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2. HIGH REDSHIFT CLUSTERS DETECTED FROM DEEP FIELDS 

To identify high redshift clusters, one has to discriminate 
the cluster member galaxies from field galaxies in deep fields. 
Note that the projection effect is very ser ious in photometric 
data. The "color cut" (iGoto et al.l 12008) or "red sequence" 
(iGladders & Yed l2005l) methods have been previously used 
for de-projection. Galaxy clusters can be detected in a certain 
redshift range when proper filters covering the4000A break 
feature are used for photometry. When the break feature shifts 
out from the chosen filters, photometric data of other bands 
must be used for color cuts. Cluster richne ss determined by 
different color cuts may not match well (e.g. lHao e t al. 2010). 
The enhanced star formation in galaxies of high redshift clus- 
ters (i.e. the Butcher-Oemler effect) also makes the cluster 
detection more difficult by the color cuts. 

The photometric redshifts of galaxies (photo-z) are deter- 
mined by comparison of the multi-band photometric data 
with the spectral energy distribution of galaxies including 
the 4OOOA break feature (e.g. ICsabai et al.l l2003l) . The 
photometric redshifts have successfully been u sed for the 
de-projection (e g., [Ei senhardt et al. 2008; Wen et al.1 120091: 
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TABLE 1 

Details of the four deep surveys used for identification of clusters. 



name 


field 


mag. 


No. of 


0-Az/(l + Z) 


photo-z 




No. of 


No. of 


No. of high-2 


No. of high-2 




(deg2) 


limit 


bands 




ref. 


this paper 


clusters 


new clusters 


clusters 


new clusters 


CFHT Wide 


35.0 


i' < 22.5 


5 


0.038 


1 


Mi, < -21.5 


631 


285 


24 


19 


CFHT Deep 


3.2 


i' < 24 


5 


0.029 


1,2 


M,, < -21.5 


202 


55 


53 


29 


COSMOS 


2.0 


i+ < 25 


30 


0.012 


3 


Mv < -20.5 


187 


127 


43 


38 


SWIRE 


33.0 


r < 24 


5-12 


0.035 


4 


Mb < -20.5 


737 


674 


116 


106 



Note. — References for photo-z: (ll lCoupon eTarU^OOgh : r21 IIlbert et mj200^ : r31 IIlbert et aPflOOgl) : (4). IRowan-Robinson et aTU20O8ft . Here, 
high-2 means 2 > 1. 



1201 ll) . When photometric redshifts of galaxies are used for 
identification of clusters, one does not have to consider the 
color of galaxies and the redshif t range. Here, we follow and 
modify the method described in IWen et al\ (l2009l) to identify 
high redshift clusters from the CFHT Wide field, the CFHT 
Deep field, the COSMOS field and the Sptizer SWIRE field. 

2.1. Photometric redshift data for deep fields 

The photometric redshifts of galaxies have been esti- 
mated for many multi-color deep su rveys, e.g., t he GOODS 
South ern Field (iMobasher et al.ll2004i) . the RCS (iHsieh et all 
I2005h . the Hubble Ultra Deep Fi eld dCoe et al.^ |2006[) . 
the Spitze r IRAC Shallow Survey (Brodwin et al. 2006|, 
CFHT (Ilbe rt et_al. |r2006), SWIRE (Rowan-Robinson et al] 
l2008h . COSMOS (Ilbert et al. 2009L the ESO Distant Clus- 
ter Survey field (felloetal.t .2009h . the AKARI survey 
(iNegrello et al.ll2009t), the Mu lti-wavelength Survey by Yale- 
Chile fcardam one etal.ll2010h . These surveys provide accu- 
rate photometries of 4 to 32 bands from UV to mid-infrared, 
so that the estimated photometric redshifts have an uncer- 
tainty of 0.01 - 0.1. For cluster detection, we restrict the 
photo-z data with a small uncertainty of <taz < 0.04(1 + z). 
In this paper, we work on the deep fields larger than 1 deg^. 
We find that photo-z data for the CFHT Wide field, the CFHT 
Deep field, the COSMOS field and the Sptizer SWIRE field 
are available for cluster identification at high redshifts (see 
details in Table [U. 

The CFHT survey is carried out with the 3.6 m Canada- 
France Hawaii Telescope in the u', g' , r' , i' and z' bands for 
three imaging sub-surveys: the CFHT Deep, the CFHT Wide 
and the CFHT Very Wide. The photo-z data of the Deep and 
Wide fields have been published and used in previous work 
(see references in Sect. 1.1). The CFHT Wide survey covers 
35 deg2 in three deep fields (Wl, W3, W4). The CFHT Deep 
survey covers 3.2 deg^ in four deep fields (Dl, D2, D3, D4). 

The COSMOS is made for a field of 2 deg^ by many ob- 
servations from X-ray to radio (XMM, Galaxy Evolution Ex- 
plorer, the Hubble Space Telescope, CFHT, United Kindon 
Infrared Telescope, Subaru, Spitzer and VLA). 

The Spitzer SW IRE is to provide photometry 
(iLonsdale et all 12003.) in the IRAC bands of 3.6 ^m, 
4.5 /Ltm, 5.8 /zm and 8.0 /.tm, the MIPS bands of 24 ^m, 70 
/im and 160 ^m for six regions: ELAIS-Nl, ELAIS-N2, 
Lockman Hole, XMM, ELAIS-Sl and Chandra Deep Field 
South (CDFS). The photometries in optical bands u, g, r, i 
and z (only here is the band, not redshift) ar e available from 
the SW IR E photometry programme and from lMcMahon et"an 
(1200 Ih and lPierreet all (120071) . 

Detailed parameters for the four fields are given in Table [T] 
including the field size, magnitude limit, number of bands for 



photometry, and the accuracy of photo-z. 

2.2. Finding high redshift clusters 

We identify galaxy clusters using the photo-z data for dis- 
crimination of lumin ous member galax ies from the four fields 
with following steps (I Wen et al.ll2009h : 

1 . Assuming that each galaxy at a given photometric red- 
shift, z, is the central galaxy of a cluster candidate, we count 
the number of luminous member galaxies of M < Afii„i, 
A^(0.5), within a radius of 0.5 Mpc and a photo-z gap of 
z ± lAAaAz- See the i\/iim and a^z values in Table [T] 
Here, the factor of 1.44 for the photo-z gap is chosen so 
that more than 85% of member galaxies of a cluster can be 
included assuming that the uncertainties of photometric red- 
shifts of member galaxies follow a Gaussian distribution. A 
larger z gap can help to include more member galaxies, but 
more background galaxies are involved. We set the abso- 
lute magnitude limit, not only to make our cluster detection 
consistent even up to high redshifts, but also to exclude less 
luminous galaxies as member galaxies which have a larger 
photo-z uncertainty. The absolute magnitudes of galaxies at 
the i'-band for the CFHT Wide and the CFHT Deep data, 
the V^-band for the COSMOS data and the i?-band for the 
SWIRE data are provided togeth er w ith pho to-z data of galaX' 
ies in their catalogs (Coupon et al.l l2009t lllbert et 



or galax- 
an i2009t 



iRowan-Robinson et aL .2008). The upper limits of absolute 
magnitude for the four fields are also listed in Table[T] The ab- 
solute magnitudes of galaxies are not available for the CFHT 
Deep field in Jlbert et al. (2006), and we estimate the ab- 
solute magnitudes in the i'-band using the K-correction of 
iFukugita et alJ 41995) according to the types of galaxies given 
bv lllbert et aiT (l2"006). 

2. We apply the friend-of-friend algorithm to the luminous 
galaxies using a linking length of 1 Mpc in the transverse di- 
rection and a photo-z gap of z ± Sctaz- We then find the 
linked galaxy with the maximum A^(0.5) as the temporary 
center of a cluster candidate. If two or more galaxies show 
the same maximum number count, the brightest one is taken 
as the temporary central cluster galaxy. The linking length 
and photo-z gap are large enough to link all member galaxies 
in the transverse distance and redshift space. This step can 
ensure to avoid multiple detection of a cluster 

3. For each galaxy with maximum A^(0.5) at z, all galaxies 
within a radius of 1 Mpc from the temporary central galaxy 
and the photo-z gap between z ± 1.44 ctax are assumed to 
be the member galaxies. The cluster redshift Zp is then de- 
fined to be the median value of the photometric redshifts of the 
recognized "members". The absolute magnitudes of member 
galaxies are re-calculated with this cluster redshift. 

4. The galaxy overdensity of a cluster, D — (A^(G.5) — 
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TABLE 2 

Clusters identified from the four deep fields. 



Name 




RA.BCG 


Decl.BCG 


Zp 


^s,BCG 


magBCG 


D 


Ngal 


R 


Ltot 


Prev. 






(deg) 


(deg) 














catalog 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


CFHT-W J02 1408.5- 


-053410 


33.53388 


-5.59269 


0.4128 


-1.0000 


18.70 


5.77 


20 


16.00 


106.41 


1,2 


CFHT-WJ021527.6- 


-053338 


33.86837 


-5.55166 


0.2848 


-1.0000 


17.70 


9.84 


25 


23.57 


120.67 


1,2 


CFHT-WJ021447.5- 


-053309 


33.69700 


-5.55293 


0.8629 


-1.0000 


20.83 


5.25 


25 


21.00 


186.48 




CFHT-WJ021242.9- 


-053219 


33.17939 


-5.54365 


1.1462 


-1.0000 


22.23 


4.55 


14 


7.43 


95.94 




CFHT-D J022435.6- 


-045505 


36.13370 


-4.91589 


0.9846 


-1.0000 


21.26 


4.30 


29 


19.76 


199.19 




CFHT-D J022723.5- 


-045424 


36.82425 


-4.91559 


0.3135 


-1.0000 


18.95 


4.28 


16 


14.31 


73.06 


1,2,3 


CFHT-D J022624.6- 


-045652 


36.59565 


-4.95964 


0.9416 


-1.0000 


21.57 


5.14 


24 


19.71 


106.97 


4 


CFHT-D J022425.3- 


-045229 


36.09802 


-4.87084 


0.9458 


-1.0000 


21.13 


4.31 


32 


21.18 


190.84 


2,3,5 


COSMOS J100313. 1+013611 


150.81335 


1.60198 


0.5119 


-1.0000 


20.25 


9.33 


15 


13.84 


49.84 




COSMOS JlOOl 17.0+013618 


150.34543 


1.62183 


0.2283 


-1.0000 


18.74 


7.25 


13 


11.37 


33.54 




COSMOS JlOOl 12.4- 


f013401 


150.28815 


1.55581 


0.3610 


0.3638 


19.08 


6.47 


36 


32.15 


105.62 


6 


COSMOS J100157.6- 


f020343 


150.49001 


2.06934 


0.4385 


0.4409 


18.93 


7.60 


14 


12.29 


44.30 


3,6 


SWIRE J003847.1- 


-433227 


9.67174 


-43.54514 


0.3609 


-1.0000 


18.88 


4.46 


14 


8.00 


28.39 




SWIRE J003847.2- 


-434822 


9.68008 


-43.79479 


0.4804 


-1.0000 


19.17 


7.57 


26 


17.11 


142.17 




SWIRE J003442.1- 


-430746 


8.67131 


-43.14276 


0.8072 


-1.0000 


21.87 


8.44 


22 


16.00 


102.18 




SWIRE J022717.1- 


-044557 


36.81958 


-4.74909 


1.3001 


-1.0000 


24.04 


4.43 


16 


11.60 


69.21 


7 



Note. — Same clusters found in different fields are repeatedly listed here for completeness with the same BCG coordinates. Column (1): 
Cluster name given by field name and J2000 coordinates of cluster center; Column (2): R.A. (J2000) of BCG; Column (3): Decl. (J2000) of BCG; 
Column (4): cluster redshift estimated from the median photo-2 of member galaxies; Column (5): spectroscopic redshift of the BCG if available. 
'-1.0000' stands for not available; Column (6): BCG magnitude at j'-band (7629 A) for the CFHT-W and CFHT-D fields, j+-band (7629 A) for 
the COSMOS field, and r-band (6230 A) for the SWIRE field; Column (7): overdensity level; Column (8): number of member galaxy candidates 
within a radius of 1 Mpc and the redshift gap of z ± 1.44 ct^^- See (taz in Table[T] Column (9): cluster richness; Column (10): total luminosity 
of cluster member galaxies, after the local background is subtracted. It is in the j'-band for the CFH T-W and CFHT- D fields, y -band for the 
COSMOS field, g-band for the SWI RE field; Co lumn (11): previous catalog containing the cluster: (1). 'Thanjavur et al.N200^ : a). IAdami eTS 
(2010); (3). Olsenet al. (2007); (4). IGrove et al. (2009); (5). Milkeraitis et al. (2010); (6). IFinogueiio v et al. (200^3) for X-ray; (7). IPierre et~ 
>2007i) for X-ray. 

This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding the form 
and content. 



limit. Here, (A^(0.5)) is the mean galaxy count of the same 
magnitude limit within 0.5 Mpc for each field for the given 
redshift gap, (Tjv(o.5) is the root mean square (rms) of the num- 
ber count (see Figure [TJ. In practices, 300 random positions 
(R.A., Decl.) are selected in the real data, and the number of 
luminous galaxies of M < A/ii,„ is counted within a radius of 
0.5 Mpc and a redshift gap between z ± 1.44 ctaz- The mean 
and the rms are statistically obtained from these 300 values. 
A larger value of D means a higher likelihood of a true cluster 
(see Table |2]l. We set the threshold, > 4, for cluster iden- 
tification. Meanwhile, we also require A^(0.5) > 8 to avoid 
the detection of poor clusters or false detection in case that the 
mean number and rms of background are very small. Now, if 
all conditions are satisfied, a cluster is identified. 

5. The real center of a detected cluster is then defined as 
the average position of luminous member galaxy candidates 
of M < Miim within radius of 1 Mpc (not 0.5 Mpc) from the 
temporary center and the redshift gap of z ± 1.44 (taz- This 
position is given in the cluster name. The cluster richness, 
R, is defined as the real number of cluster galaxies, which 
is obtained from the number of member galaxy candidates, 
A'gai, within a radius of 1 Mpc from the real center of a cluster 
and the redshift gap subtracted by the local contaminations 
(A'cb) of foreground and background galaxies, so that R = 
Ngal — (Nch)- The contamination, (Nch), is estimated from 
the local distribution of th e galaxies in this redshift gap (see 
details in lWen et al.ll2009l) . 

6. The total luminosity of each cluster, itot, is calculated 
as the total luminosity of member galaxies within the region 
after the similar subtraction of the average background con- 
tribution. 
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radius of 0.5 Mpc and a redshift gap between z ± 1.44 ctaz ■ 
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(A^(0.5)))/(Tjv(o.5)^ is estimated from the member number 
and the background galaxy density of the same magnitude 
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Fig. 2. — Examples for clusters in the CFHT Deep field ident ified only 
from d ata of Ilbert et al. (2006) (upper) and only from data of Cou pon et alj 
(200P) (bottom). The images have a field of view of 1.5' X 1.5'. The most 
of bright red galaxies in the images are member galaxies. A color version of 
this figure is available in the online journal. 
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Fig. 3. — Redshift distribution of detected clusters from the four deep 
fields. The grey area indicates the previous known clusters in our sample. 



As listed in Table |2] we get 63 1 clusters from the CFHT 
Wide field, 187 clusters from the COSMOS field and 737 
clusters from the SWIRE field. In t he CFHT Dee p field 
two sets of photomet ric redshifts by ICoupon et al.l (l2009b 
and lllbert etal] (|2006|) are available for galaxies with simi- 
lar uncertainties. Both are used for our cluster detection in- 
dependently. We get 163 clusters using the photo-z data of 
ICoupon et a l. (2009) and 105 clusters using the photo-z data 
of fllbert eTaiJ (12006). Among t hem, 66 clusters are d etected 
from both, 97 clusters only fr om [Cou pon et al1(l2009b and 39 
only from lllbert et aP (l2006l) . Figure [2] shows tvyp ex ample 
clusters identified only from the data of Ilbert et al. ('2006) and 
two examples only from the data of Coupon et al. (2009). All 
clusters have an overdensity D > A m one dataset but less 
than 4 in the other Merging two samples gives 202 unique 
clusters from the CFHT Deep field. 

We notice that there are the overlapped regions between the 
four fields: the COSMOS field is overlapped with a part (D2) 
of the CFHT Deep field; a part (Wl) of the CFHT Wide field 
is overlapped with the Dl region of the CFHT Deep field and 
the XMM region of the SWIRE field. Therefore, some clus- 
ters are repeatedly detected from the photometric redshift data 
of different fields. Merging 1757 entries of clusters in Ta- 
ble |2] from the four fields gives 1644 clusters in total. Among 
them, 1088 clusters are newly identified, and more than half 
of them are identified from the large SWIRE field which has 
not been searched for clusters previously. Among 228 clus- 
ters of z > 1, 191 are newly identified, and again, more than 
half of them are identified from the SWIRE field. See Tabled 
for the numbers of clusters identified from each field. 

Figure |3] shows the redshift distribution of detected clus- 
ters in each field in the range of 0.1 < z < 1.6. The me- 
dian redshifts of the clusters from the CFHT Wide, the CFHT 
Deep, the COSMOS and the SWIRE are 0.63, 0.86, 0.85, 



0.70, respectively. According to the redshift distribution of 
clusters, clusters of z < 0.7 (i.e. less than the peak redshift) 
in the CFHT Wide field probably have their luminous member 
galaxies selected with a similar completeness. At higher red- 
shifts, the richness is biased to a smaller value, since only the 
brighter galaxies can be detected because of the observational 
magnitude limit. Similarly, the redshift is z < 0.9 for clus- 
ters in the CFHT Deep field, z < 1 for the COSMOS field, 
z < 0.8 for the SWIRE field. 

We check the companions within a radius of 2 Mpc and 
z ± 3 (JAz of the detected clusters, and find that 76% of clus- 
ters have no companion, 20% of clusters have one companion, 
and 4% of clusters have two companions. Companions for a 
small number of clusters probably indicate superclusters or 
the cluster mergers. 

2.3. False detection rate 

Because of the projection effect and the uncertainty of 
photo-z, it is possible that a false cluster is detected by chance 
in the photo-z data of galaxies by the procedures above. We 
have to estimate the false detection rate of our cluster detec- 
tion algorithm. 

Following the method of IWen et akl (l2009h . we perform 
Monte Carlo simulations based on the real data for the pur- 
pose. First, each galaxy in the real dataset is forced to have 
a random walk in the sky plane towards a random direction 
with a step length of a random value larger than 1.0 Mpc but 
less than 2.5 Mpc. Second, we shuffle the photo-z and ab- 
solute magnitude of all galaxies in each field. These proce- 
dures should eliminate real clusters but reserve the statisti- 
cal properties of original dataset. We apply the cluster detec- 
tion algorithm to such a 3-D shuffled data, to see how many 
"false clusters" can be detected with the criteria of £> > 4 
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Fig. 4. — Distribution of tlie overdensity of detected clusters from real data 
and the "false" clusters from the 3-D shuffled data (grep area). 




Spectroscopic redshift, 

Fig . 5 . — Comparison between the cluster redshift Zp (the median of photo- 
z of member galaxies) and the spectroscopic redshift of BCG, z^: dots for 48 
CFHT clusters with a standard deviation, o"(^^_^^jy(^^^^), of 0.017, open 
circles for 48 COSMOS clusters with a standard deviation of 0.006, and stars 
for 26 SWIRE clusters with a standard deviation of 0.018. 



and A^(0.5) > 8. We repeat such tests 100 times, and get 
the false detection rate of 9.9%, 2.9%, 3.2% and 7.9%, re- 
spectively, for the CFHT Wide data, the CFHT Deep data, the 
COSMOS data and the SWIRE data. Figure|4]shows the aver- 
age overdensity distribution of the the detected clusters from 
the four fields, together with that of "false clusters" from the 
3-D shuffled data. The "false" clusters preferably have a low 
overdensity. 

2.4. The brightest cluster galaxies 

We recognize the BCG of a cluster as the brightest galaxy 
within a radius of 0.5 Mpc from the cluster center and a photo- 
z gap of z ± 1.44 (Jaz- Furthermore, we check the images of 
all clusters and BCGs to clean out any possible contamination, 
e.g., blue spiral galaxies mis-identified as BCGs and bad pho- 
tometries. For the CFHT Wide and Deep clusters, we inspect 
the color composite CFHT images at the website'. For the 
COSMOS and SWIRE clusters overlapped in the CFHT field, 
we view the CFHT images, otherwise the COSMOS images^ 
or the SWIRE images'*. We also inspect the SDSS color com- 
posite images'* for clusters of z < 0.6 located at the SDSS 
field. The coordinates and the AB magnitudes of the BCGs 
are given in Table |2] The AB magnitudes are taken from the 
data for the CFHT, the COSMOS and the XMM region of the 
SWIRE field. The Vega magnitudes are taken from that for the 
SWIRE field (except the XMM region) and are transfer to AB 

' http://ww wl.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/community/CFHTLS-SG| 
/docs/csky.html 

2 http://irsa.ipac.caltech.edu/data/COSMOS 
^ http://irsa.ipac.caltech.edu/data/SPITZER/SWIRE 
http ://cas . sdss . org/astrodr7/en 
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Fig. 6. — Correlations between the absolute magnitudes (i.e. the luminos- 
ity) of BCGs with the cluster richnesses and the total luminosities. 



magnitude^ for the rest SWIRE field via tab = '"Vcga + 0.17. 

Because BCGs are brighter than other member galax- 
ies, some of them have the spectrosco p ic red shifts ob- 
tained alread y, e.g. Le Fevre et al.l (l2005l) for the 
CFHT field, iLinveta l JIOOTh for flie COSMOS field. 



iRowan-Robinson et alj (12008') for the SWIRE field. We also 
check the spec troscopic red shifts of BCGs from the SDSS 
DR7 database ("Abaz ajianet alJ 1200 9). Figure |5] shows the 
difference between the median photometric redshift of clus- 
ter members and the spectroscopic redshift of BCGs, [zp — 
Zs)/{y + Zs). The standard deviations roughly indicate the 
accuracy of redshift estimate for clusters. 

The luminosities of BCGs are more or less correlated to 
the cluster richness and the cluster total luminosities (see 
Fig. |6ll. Richer clusters tend to have more luminous BCGs. 
The strongest correlation appears between absolute magni- 
tude of the COSMOS BCG and total luminosity. 

2.5. Matching for different detections 
' |http://www.sdss.org/dr7/algorithms/sdssUB VRITransfomi.html I 
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Fig. 7. — Comparison between the richnesses of matched clusters in the 
overiapped regions between the CFHT Wide and the SWIRE fields and be- 
tween the CFHT Deep and the COSMOS with different absolute magnitude 
limits. 



Note that the absolute magnitude limits of the member 
galaxy selection are different for the four deep fields (see Ta- 
ble [TJ. That is Me = -21.5 for the CFHT field, My = 
-20.5 of an aperture of 3" for the COSMOS field, and 
Mb = -20.5 for the SWIRE field. Different Hmits allow 
the member galaxies of different luminosities to be selected 
in the cluster detection algorithm. The richness of a cluster 
therefore depends on the absolute magnitude limits. 

Clusters in some overlap regions have been detected from 
the dataset of two deep fields via the cluster detection al- 
gorithm, matched within a separation of 1 Mpc and a red- 
shift difference of Az < 0.04(1 + z). We compare their 
richnesses obtained from the two field datasets with differ- 
ent absolute magnitude limits at different bands, if luminous 
member galaxies are almost completely selected. We get 6 so 
matched clusters in the overlapped region of the CFHT Wide 
field and the CFHT Deep field, 24 matched clusters in the 
COSMOS and CFHT Deep fields, 31 matched clusters in the 
SWIRE and CFHT Wide fields. Figure [T] compares the rich- 
nesses of the matched clusters obtained from the overlapped 
regions between the CFHT Wide and the SWIRE fields and 
between the CFHT Deep and the COSMOS fields. As ex- 
pected, these richnesses are correlated. 

As shown in the Figure |3] and mentioned in Sect. 1.1, a 
large amount of clusters have previously been detected in the 
CFHT Wide and Deep fields. For example, Olsen et al. (2007|) 
applied a matched-filter algorithm to the i'-band four CFHT 
Deep fields, and identified 162 cluster can didates with a n esti- 
mated redshift uncertainty of ctaz = 0.1. iGrove et"ail (l2009h 
applied the method to the r' and z' bands the four CFHT Deep 
fields, and found 1 14 and^247 cluster candidates from the two 
bands, respectively. iThaniavur et al.l (120091 and private com- 
munication) used a red sequence method to the CFHT Wide 
field of the 161 deg ^ data and found 5804 cluster candidates. 
lAdami et al] (1201 Oh used a 3-D method and identified 1029 
clusters from the 28 deg^ CFHT Wide field and 171 clusters 
from the 2.5 deg ^ CFHT Deep field with a mass greater than 
1.0 X 10" Mq. iMilkeraitis etaP (12010 ) used 3D-Matched- 
Filter to identify clusters from the CFHT Deep fields, and 
found 673 clusters in the redshift range of 0.2 < z < 1.0. 
We merge these previous cluster samples, and get 760 unique 
clusters in the CFHT Deep field, and 5848 unique clusters 
from the CFHT Wide field. We match those with our sam- 
ple by a separation of 1 Mpc and a redshift difference of 
Az < 0.2, and get 346 matched clusters in the CFHT Wide 
field, and 144 matched clusters in the CFHT Deep field. In 
Figure [8] we compare the distribution of known clusters and 
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Fig. 8. — Known clusters and our detections. Upper panel: the percentage 
of our detected clusters are previously known in the CFHT Deep field (open 
circle), especiall y by 3 -D methods (black dots) in Adami et al. (2010) and 
iMilkeraitis et ai] 120101) . Lower panel: the number of the 3-D clusters we 
detect. Apparently we detect the 3-D clusters of only a high overdensity. 
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Fig. 9. — Comparison of the richnesses of matched clusters we define by 
the number of luminous galaxies with those given in Olsen et al. ( 2007) and 
'Grove et al.' ^2009) defined as the average number of L* galaxies based on a 
single band photometric data. Here L* is the characteristic luminosity in the 
Schechter luminosity function. 



our detections. At low redshifts, almost all clusters we detect 
are known previously. However, we detect some new clusters 
at high redshifts which are not previous ly known. Among 
the clusters detected b y the 3-D method (lAdami et al.ll20I0t 
IMilkeraitis et al.l20IOl) . we detect only those rich ones in gen- 
eral which have a high overdensity of D > 4. 

In Figure|9] we compare our richnesses of matched clusters 
of z < 0.9 in the C FHT Deep field with t he richness given by 
lOlsen et al.l (|2007|) and iGrove et al] (12009) as the equivalent 
number of L* galaxies, based on photometric data of a single 
band, here L* is the characteristic luminosity in the Schechter 
luminosity function. Though with a large scatter, richer clus- 
ters found by previous authors show a larger richness by our 
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Fig. 10. — Correlations of the X-ray luminosity with the cluster richness, 
^COSMOS, (the middle panel) and the total V-band luminosity, ^'qqsmoS' 
(the right panel) with for 24 COSMOS clusters in the redshift range of 0.2 < 
z < 1.0 (see the left panel). 



definition. 

2.6. X-my emission of clusters 

Some clusters in the COSMOS field ha ve X-ray emission 
detected. Using the XMM-Newton data, Finog uenov et alj 
12007) identified 72 X-ray clusters in the 1.7 deg^'of the COS- 
MOS field, 24 of which are matched with our sample within a 
projected separation of rp < 1 Mpc and a redshift difference 
of Az < 0.05 2(T of the cluster redshift uncertainty in 
iFinogu enov et al.| |2007h . These X-ray clusters have a redshift 
in the range of 0.2 < z <1 (see Figure [Toll. 

Similar to low-redshift galaxy clusters (IWen et al.l l2009l) . 
significant correlations are found between the X-ray luminos- 
ity of clusters and the cluster richness {R, the middle panel of 
FigurefToli or the total l/-band luminosity (Lv, the right panel 
of FigurefToli. The best fittings to the data give 

^ (_o.75± 0.43) + (1.62 ±0.34) log i?cosMos, 

(1) 

and 

logL^^42 ^ (_2.02± 0.48) + (1.75 ±0.26) log L^^°Mos: 

(2) 

where i^'"*^ refers to X-ray luminosity in the 0.1-2.4 keV 
band in unit of 10^^ erg 
t/-band luminosity in unit of 10 ijq 
measurements. The correlations suggest that the richness we 
define closely relates to the properties of clusters. 



^COSMOS refers to the total 
^° from the COSMOS 



3. STELLAR POPULATION AND COLOR EVOLUTION OF BCGS 

The colors of BCGs provide important information on their 
stellar population. Since the clusters in our sample have large 
redshifts, one can easily study the stellar population and color 
evolution of these BCGs. Here, we work on color evolution 
as a function of redshift, in the light of a stellar population 
synthesis model. 

For the CFHT data, the BCG magnitudes at u'g'r'i'z' bands 
are available from Coupon et al. (2009). The r' (7480 A) and 
z' (8930 A) bands correspond to the rest-frame optical wave- 
length for observed BCGs up to redshift z ~ 1, and the g' 
band (4860 A) corresponds to the rest-frame optical wave- 
length at low redshifts and the UV wavelength at high red- 
shif^ts. The BCG magnitudes at the r', z' and g' bands are 
used to define the BCG colors for the clusters in the CFHT 
Wide and Deep fields, as shown in Figure[TT]and Figure [T2l 

To explain the color changes, we show the evolution lines 
calculated by using the stellar population synthesis models 
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Fig. 1 1 . — The colors, g' — z',r' — z' as a function of redshift for BCGs in 
the CFHT Wide field. Three lines represent the color evolution with redshift 
from a model of passively evolved galaxies fonned at redshift Zf = 1 (solid 
line), Zf = 2 (dashed line) and Zj = 3 (dotted line), respectively. To outline 
the evolution at lower redshifts (z < 0.6), the same-color data of 500 BCGs 
from the SDSS clusters in iWen et al.i Ii2009i) are included in the plots. 



dBruzual & Chai-lot"2003', BC03, herea fter). We ado pt the 
stellar evolution tracks of Padova 1 994 (iGirardi et al.|[l996i) . 
the Basel3.1 stellar spectral lib rary ([Westera et al.ll2002l) . and 
the initial mass function of Chabrie772003^ and the Solar 
metallicity. The solid line in Figure [TT] stands for the color 
evolution for a galaxy formed at redshift Zf = 1 and evolved 
passively with cosmic time, and the dashed and dotted lines 
for galaxies formed at redshift z/ = 2 and 3, respectively. We 
find that the BCG color*", r' ~ z' , is consistent with the BC03 
model of Zf = 2 and 3, but significant redder than that of 
Zf = 1. The BCG color, g' — z' , of z < 0.8 is consistent with 
models of z/ ~2 and 3, but significant bluer than the models 
at higher redshift. 

For the COSMOS data, the B (4440 A), r+ (6232 A) and 
m-s.etim bands are used to define the BCG colors. The B 

l (l2009l) . 



and r+ band data are available from lllbert et al.l (120091) . The 
B band corresponds to the rest-frame optical wavelength at 
low redshifts, but the rest-frame UV wavelength at high red- 
shifts. We obtain the 3.6/im flux within an ap erture of 1.9" 
diameter from the Spitzer S-COSMOS survey (ISanders et al.l 
lIOOTh . Since the optical fluxes are measured over an aper- 
ture of 3" in diameter which encloses 75% of the flux for a 
poi nt-like source, we therefore convert the 3.6/im flux follow- 
ing Slberretan (2009) to match the aperture of optical data. 
The AB magnitude at 3.6/im is then calculated via equation, 

^ When comparing the color data with model, we have to systematically 
add O.I to r' - z' and 0.2 to g' - z' of BCGs in the CFHT Wide field 
(see Figure rm to get consistence a t low redshift (z < 0.3). For the BCGs 
from the CFHT Deep field in Figure [T2] we have to systematically add 0.2 to 
r' — z' and 0.3 to g' — z' . These calibrations are necessary to compensate 
the possible systematical bias in photometric data of different bands. 
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Fig. 12. — The same as Figure [171 but for BCGs of clusters found in the 
CFHT Deep field. 
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Fig. 13.— The color evolution of BCGs in the COSMOS field. The lines 
are calculated from the BC()3 model as the same as Figure 1111 but for the 
colors, r+ — ma.e^m and B — ms.e^im. 



23.9 - 2.5 log(/3.6^m), where /s.eMm is in unit of 

Figure [T3] shows the color evolution of BCGs in the COS- 
MOS field. We find that the color, r' — 3.6/iTO, is consistent 
with the BC03 model for Zf ~ 2 and 3. The color data of 
B — TO3.6Aim at z < 0.5 is consistent with that of models of 




Redshift, z 

Fig. 14.— The color evolution of BCGs in the Spitzer SWIRE field. The 
lines are calculated from the BC03 model as Figure [171 but for the color of 



Zf =2 and 3, but significant bluer than the models at higher 
redshift. 

For the SWIRE data, the magnitudes of r and 3.6Aty» bands 
are available from lRowan-Robinson et all (1200 8*) and are used 
for color evolution study (see Figure [Mt. The result is consis- 
tent with that from the COSMOS clusters. 

The wavelenghts of r', z' , r+ and ms g^m bands corre- 
spond the optical or infrared band in the rest-frame, and 
the observed emission of galaxies in these bands is domi- 
nated by old stellar populations. The color evolution shown 
by these high redshift data shows that the stellar population 
in m ost BCGs have been f orme d at z > 2, in agreement 
with IStott et all (e.g., l2008h and IWhilev et all (|2008). The 
data of colors g' — z' and B — m^ ^^rn at high redshifts in 
Figure [TT1 - [T3] are related to the enhancement rest-frame UV 
color which indicates star formation in the high redshift BCGs 
JPonahue et al.,,2010; Hicks et al. 2010). 

4. CONCLUSIONS 

Using the photo-z data of the four deep fields, we identified 
631 clusters, 202 clusters, 187 clusters and 737 clusters, re- 
spectively, from the CFHT Wide field, the CFHT Deep field, 
the COSMOS field and the SWIRE fields. Clusters are rec- 
ognized when an overdensity is Z? > 4 and more than eight 
luminous galaxies are found within a radius of 0.5 Mpc and 
a photo-z gap. These clusters have redshifts in the range of 
0-1 z < 1.6. Because of overlapping areas among the 
four fields, merging these cluster samples gives 1644 clusters, 
of which 1088 clusters are newly identified and 228 clusters 
have a redshift z > 1.0. The false detection rate is estimated 
to be less than 10%. The cluster redshift is estimated as the 
median photo-z of member galaxies. Richer clusters tend to 
have more luminous BCGs. The cluster richness and total lu- 
minosity are tightly related to the X-ray luminosity. 

The BCG colors and their evolution are studied for the large 
sample clusters in the four fields. We compared them with a 
passive galaxy evolution model using the stellar population 
synthesis. The color evolution is consistent with the BC03 
model in which the stars in BCGs are formed at z/ > 2 and 
evolved passively. The systematical enhancement rest-frame 
UV color indicates star formation in these BCGs. 
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